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We derive an exact general axi-symmetric solution of the coupled gravitational and elec- 
tromagnetic fields in the tetrad theory of gravitation. The solution is characterized by four 
parameters M (mass), Q (charge), a (rotation) and L (NUT). We then, calculate the total 
exterior energy using the energy-momentum complex given by M0ller in the framework of 
Weitzenbock geometry. We show that the energy contained in a sphere is shared by its 
interior as well as exterior. We also calculate the components of the spatial momentum to 
evaluate the angular momentum distribution. We show that the only non-vanishing compo- 
nents of the angular momentum is in the Z direction. 



1. Introduction 



O . 

A general stationary axi- symmetric object, with gravitomagnetic monopole and dipole 
• • ! moments associated with nonzero values of the NUT and Kerr parameters L and a respec- 
. ^ | tively, is described by Kerr-NUT (Newman- Unti-Tamburino) spacetime [1] which is a useful 
^ ; model for exploring gravitomagnetism [2] . The Kerr-NUT spacetime and its spacial cases are 
c5 'all belong to larger class of stationary axi-symmetric type D. Carter [3] has found vacuum 
solutions of Einstein equations for which the Hamilton- Jacobi equation of the geodesic is 
separable. The stability of the general axi-symmetric spacetime is probed by studying their 
perturbation by fields of various spin [4, 5]. 

Einstein's general relativity (GR) is a very successful theory in describing long distance 
phenomena. However, it encounters from serious difficulties on short distances. The main 
problem appears in all attempts is the quantization [6, 7]. Also, the Lagrangian structure of 
GR differs from the ordinary microscopic gauge theories. In particular, a covariant conserved 
energy-momentum tensor for the gravitational field can not be constructed in the framework 
of GR. Consequently, the study of alternative models of gravity is justified from the physical 
as well as from the mathematical point of view. Even in the case when GR is a unique true 
theory of gravity, consideration of close alternative models can shed light on the properties 
of GR itself. 

It is well known that in GR an energy-momentum tensor fails to satisfy some certain 
conditions [8]. This is usually related to the equivalence principle which implies that the 
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gravitational field can not be detected at a point as a covariant object. This can be viewed as 
a purely differential-geometric fact. Since the components of the metric tensor are managed 
by a system of second order partial differential equations therefore, the energy momentum 
quantity has to be a local tensor constructed out from the metric components and their 
first order derivatives. The corresponding theorem of (pseudo) Riemannian geometry states 
that every expression of such type is trivial. It is natural to expect that this objection for 
the existence of a gravitational energy-momentum tensor is directly related to the geometric 
properties of the (pseudo) Riemannian manifold. This objection can be lifted in alternative 
model, even connected with the geometry of the manifold. 

In recent time teleparallel structures on spacetime have evoked a considerable interest for 
various reasons. They were considered as an essential part of generalized non-Riemannian 
theories such as the Poincare gauge theory ( [9] ~ [13] and references therein) or metric-affine 
gravity [14]. Physics relevant to geometry may be related to the teleparallel description of 
gravity [15, 16]. Teleparallel approach is used for positive-gravitational-energy proof [17]. 
The relation between spinor Lagrangian and teleparallel theory is established [18]. 

M0ller has shown that the problem of the energy-momentum complex has no solution in 
the framework of gravitational field theories based on Riemannian spacetime [19]. In a series 
of papers, [19, 20, 21] he was able to obtain a general expression for a satisfactory energy- 
momentum complex in the absolute parallelism space. The Lagrangian formulation of the 
theory was given by Pellegrini and Plebanski [22]. Quite independently Hayashi and Nakano 
[23] formulated the tetrad theory of gravitation as a gauge theory of the spacetime translation 
group. In these attempts, the admissible Lagrangians were limited by the assumption that 
the field equations has the Schwarzschild solution. M0ller [24] abandoning this assumption 
and look for a wider class of Lagrangians by constructing a new field theory. His aim was to 
get a theory free from singularity while retaining the main merits of GR as far as possible. 
Meyer [25] showed that Miller's theory is a special case of the Poincare gauge theory [26, 27]. 
Sdez [28] generalized M0ller theory into a scalar tetradic theory of gravitation. 

The tetrad theory of gravitation based on the geometry of absolute parallelism [22]~[34] 
can be considered as the closest alternative to GR and it has a number of attractive features 
both from the geometrical and physical viewpoints. Absolute parallelism is naturally for- 
mulated by gauging spacetime translations. Translations are closely related to the group of 
general coordinate transformations which underlies GR. Therefore, the energy-momentum 
tensor represents the matter source in the field equation for the gravitational field is just like 
in GR. 

The tetrad formulation of gravitation was considered by M0ller in connection with at- 
tempts to define the energy of gravitational field [24, 35]. For a satisfactory description of 
the total energy of an isolated system it is necessary that the energy-density of the gravita- 
tional field is given in terms of first- and/or second-order derivatives of the gravitational field 
variables. It is well-known that there exists no covariant, nontrivial expression constructed 
out of the metric tensor. However, covariant expressions that contain a quadratic form of 
first-order derivatives of the tetrad field are feasible. Thus it is legitimate to conjecture that 
the difficulties regarding the problem of defining the gravitational energy- momentum are 
related to the geometrical description of the gravitational field rather than are an intrinsic 
drawback of the theory [36, 37]. 

The main aim of the present paper is to derive a general axi-symmetric solution in 
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the tetrad theory of gravitation for the coupled gravitational and electromagnetic fields. 
Using this solution we then, calculate the energy and angular momentum using the energy- 
momentum complex given by M0ller [24] and Mikhail et.al [38]. In §2 we derive the field 
equations of the coupled gravitational and electromagnetic fields. In §3 we obtain a new 
exact analytic axi-symmetric solution characterized by four parameters in the tetrad theory 
of gravitation. In §4 we calculate the energy and angular momentum distribution of this 
solution using M0ller's energy-momentum complex [24]. The final section is devoted to 
discussion and conclusion. 



2. Tetrad theory of gravitation and electromagnet ism 



In the Weitzenbock geometry the fundamental field variables describing gravity are a 
quadruplet of parallel vector fields [26] e/*, which we call the tetrad field in this paper, 
characterized by 

D v e i » = d v e i ' i + r' i Xv e i x = 0, (1) 

where T^\ u defines the nonsymmetric afline connection coefficients. The metric tensor g^ u is 
given by 

with the Minkowski metric rjij = diag(+l , — 1 , — 1 , —1) L The curvature tensor defined by 
T x fiu , given by Eq. (1) i.e. R a ^(T) , is identically vanishing [26]. 

The gravitational Lagrangian L G is an invariant constructed from g^ and the contorsion 
tensor 7^ given by 

li^up 6 ^ejj, ; p , (2) 

where the semicolon denotes covariant differentiation with respect to Christoffel symbols. 
The most general gravitational Lagrangian density invariant under parity operation is given 
by the form [23, 26, 38] 

C(e) G = bL G (e) = e (a^, + a 2l ^ p llivp + a 3 7^ P 7^) (3) 
and $ M being the basic vector field defined by 

*m = 7V (4) 

where e = det(e a p). Here ol\ , a 2 , and 0:3 are constants determined such that the theory 
coincides with general relativity in the weak fields [23, 24]: 

ai = --, a 2 = - , a 3 = -(1 - A) , (5) 

K K K 



*Latin indices (i,j,k,- ■ •) designate the vector number, which runs from (0) to (3), while Greek indices 
(/x, v, p,- ■ •) designate the world-vector components running from to 3. The spatial part of Latin indices is 
denoted by (a, b, c, ■ ■ •), while that of Greek indices by (a, (3, 7, • • •). 

t Latin indices are rasing and lowering with the aid of rjij and rfK 
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where k is the Einstein constant and A is a free dimensionless parameter-^. The vanishing of 
this parameter makes the theory coincides with GR formulated in teleparallel geometry. 

The electromagnetic Lagrangian density L e . m , is [32] 

L e .m. = -^ p g UCJ F^F p(T , with F^ being given by = d^A v - d u A^ (6) 

where A^ is the vector potential. The variation of the gravitational Lagrangian given by Eq. 
(3) with respect to the tetrad field gives one equation of motion. M0ller assumed that the 
energy- momentum tensor of matter fields, i.e., T^ v is symmetric. Therefore M0ller divided 
the resulting field equation into symmetric part and skew symmetric part [24]. 

The gravitational and electromagnetic field equations for the system described by Lq + 
L e m . are the following: 

G^v + H^ v = —kT^v , (7) 
Kp, = , (8) 

d v (V^gF^) = (9) 
with G^y being the Einstein tensor of GR. Here H^ v and K^ v are defined by 



jj del. , 
tl^v = A 



(10) 



and 

= A $^ jV — <& u )/U — $p (^f P /iu ~~ ^ P uji) + l^v p . p ) (11) 

and they are symmetric and antisymmetric tensors, respectively. The energy-momentum 
tensor T^ v is given by 

Fjxv = -g p(J F w F va + - A g, u F p °F pa . (12) 

It can be shown [26] that the tensors, H^ v and K^, consist of only those terms which 
are linear or quadratic in the axial- vector part of the torsion tensor, a M , defined by 

a„ =■ \e^l up % (13) 

where e MJ , pcr is defined by 

tfivpo y/ g&fivpcj i (14) 

with S^pa being completely antisymmetric and normalized as <5 i23 = — 1- Therefore, both 
and F^ u vanish if the a M is vanishing. In other words, when the a M is found to vanish 
from the antisymmetric part of the field equations, (8), the symmetric part (7) coincides 
with the Einstein equation formulated in the Weitzenbock spacetime. 



-^Throughout this paper we use the relativistic units , c = G = 1 and k = 8ir. 
§ Hcaviside-Lorentz rationalized units will be used. 
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3. Exact Analytic Solutions 



In this section we will seek for a solution satisfying the following conditions: The parallel 
vector fields having the form 



6* + Mld k + 



k . Q 2 - L 



2 M 



(15) 



Here M, Q 2 and L 2 are free parameters, / M and m M are quantities satisfying the conditions 



■q^m^my = 0, 



fU = 0, 



(16) 



and l k and m fc are defined by 



k def. rfc «i/ 



(17) 



Applying (15) to the field equations (7)~(9) we obtain the values of 1^ and in the form 



to = v ^ 



T + r 2 + /i 2 



Q I 



— ^ e a/ 33aa; p 



m. 



(18) 



where 



n 

1 



T = J(r 2 -a 2 ) 2 + 4a 2 z 2 , Q 



I r 2 — a 2 + T 



(19) 



with r = v^ 2 + y 2 + z 2 and e a/ g 7 are the three dimensional totally antisymmetric tensor with 
^123 — 1 and a is the angular momentum per unit mass. The form of the electromagnetic 
potential has the expression 



A u = — - VE L, and the Maxwell field F av — 



q 



■ (20) 



where S and 1^ are defined in Eqs. (18) and (19). 

For the solution given by Eqs. (18) and (20), the axial vector part a^ of the torsion 
tensor vanishes, 

a„ = 0, (21) 
and the metric is identical to Kerr-NUT metric in GR. 

Writing explicitly the tetrad (15) using (18) and (19) we obtain* 

(2Mp - Q 2 + L 2 )p 2 



3 (0) o = 1 



2pi 



*We will denote the symmetric part by ( ), for example , A^ v ) = (l/2)(A tlI/ +A vll ) and the antisymmetric 
part by the square bracket [ ], A^ u] = (1/2)(A^ - A Vfl ) . 
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3 (0) 



, a ... a 2 z x 3 1 2Mp-Q 2 + L 2 p 4 m 
P P 2 P 2 Pl (p 2 + a 2 ) 



{Ob CL^ ( CI 

x l x p - 2-e k3( pX l) x k + — e k l3 e m p 3 x k x m + zl {px l - -e k3 l x k }5 p 

\ (2Mp -Q 2 + L 2 )p A 



+{pxp - -^e kW x k }5 13 



V^^ 3! | 2 Pl (p 2 + a 2 ) 2 



(22) 



with S l p being the Kronecker delta [32] and 



pi = p 4 + a 2 z 2 . 



Solution (22) with Eq. (20) satisfy the field equations (7)~(9) and the associated metric has 
the following form 



(2Mn -a 2 -Q 2 + L 2 ) ri 2 
r*! 4 + a?z 2 



ds 2 = dt 2 -dx 2 -dy 2 -dz 2 



zdz r\(xdx + ydy) + a(xdy — ydx) 



dt+ + 



ri 2 + a 2 



(23) 



with r\ is the radial parameter related to the Cartesian radius r by 

4 / 2 2\ 2 22 n 

ri — (r — a )ri — a z = 0. 



4. The energy and angular momentum associated with the 
general axi-symmetric solution 



The superpotential is given by [24, 38] 
(-a) 1 / 2 

W/ A = K -^—P XP J vX Wf x 9^ ~ \9r»l xp ° - (1 - 2\)g T ^ PX ] , (24) 



where P XP J vX is 



= 8/9^ + 5 p T g a / x - 5/g xp vX (25) 



with g pa uX being a tensor defined by 

9po vX d = S p »5 a x - 5/5 p \ (26) 
The energy-momentum density is defined by [24] 

V=ZV\ A . (27) 



The energy E contained in a sphere with radius R is expressed by the volume integral [39] 



E(R)=J r=R J jr d 3 x. (28) 

For convenience of the calculations, we work for small values of the rotation parameter a, 
so we will neglect the terms beyond its fourth order. With this approximation, we have the 
covariant components of the parallel vector fields (22) as 



a 4 \( 4(Q 2 -L 2 )\ ( 24(Q 2 

- 3M - - 30M ^- 

8r 5 y r J y r 



r J \ r 



A(Q 2 -L 2 )\ ( 2A(Q 2 -L 2 )\z 2 ( 24(Q 2 - L 2 
3M - '- - 30M y — '- — + 35M ^ = 

r / q 2 -l 2 



/ V 

r r 2 y 2r J a/3 



2 )\ z A 



I r / r 2 a/33 8r 5 | 2r 



I r a ' 2r 4 \ r I 

,(g!^! + f 8M - 7 W 2 - L2 ')4- 

j 2r y r y r 2 

^-fi2M- 7 ^ 2 ;^UU 



21(Q 2 



2r 



r 

L 2 ) \ : 



3 Ll v 2r 

' r2\ 



^|n Q -2|2^2M-^ 



m n /"V 2a fl \ / r Q 2 - L 2 \ a 2 [l f / 3(Q 2 -L 2 )U 2 1 , 

e( °« = 5 '« + {-f - ^W') [ M - J " ^ [2 \ M + y 5M " r j H J ^ n ° 

" (" " ^) - (m - S^) + £ [{. (m - S^) £ 

^ - (m - S^j ^ w] - £ [{ (- - ^ 



ivii I ji" 



^ / 40(Q 2 -L 2 
14M— + 63M ^ 



_ f 7M - iigl^l) 4}^- - J (M - - f 8M - 

y r / r J r [ y r /\ r 

+ u - 4k + J f m - - u - 4 

y r y r 4 J [ y r /\ r I r 

and the contravariant components of the parallel vector fields (22) have the form 



e (0 ) = 



1 / , Q 2 -L A \ a" 
1 + - M- + — 

r 2r / 2r 3 



M -g^!)-( 3 M- 2 W 2 - L2 ))4 

r ) \ t It 



+ 



a 

8r5 



3M 



4(Q 2 



30M 



24(Q 2 



e (0) Q = S al b^ a = -b (a) °, 

en f n l na 2 « 



v 7 ' 



a 2 



~2 



a* 

r 4 



^ + } - w» - (m - 9^ ?V + £ [{ (5M - 



y H y r 
2^ ral3 *,r Q 2 — L 2 \ 



n -{ M - J T n ^'H + 8^ [\ [ 5M ~ 

y y ryr 2 y r y r 4 



£ 2 ) U 2 



r 



Also with this approximation the covariant components of the metric tensor have the form 

r9\ 9 T/ ^9 7- 9 \ / n / ^~>9 r9\\ 91 



(30) 



5oo = 1 - - M 7, - -3 \ M 



a* 



90a — — 



'ITT ° 1V1 

-{m-1 



2r 

( A(Q 2 - L 2 
3M- — 



lenxs 01 xne mexric xensor nave 

-)-(sm-HM)£!] 



2r 



L 2 \ 2n„ 2o 



+ y 2 ))n a -(M-9LJl 



( 2A(Q 2 -L 2 )\z 2 ( 2A(Q 2 -L 2 

30M K — — + 35M ^ 

\ r J r \ r 

\ 2r 



2 )\ z 4 

r 4 



/? 2a 2 
e a /33« p + — 



1 2 1 2 



Z A 3 

r 



4r 5 L [ 2r \ r ) r 2 \ 

_ 2 { 2 ( 2M _w^!))_( 12M _ZW!_i! 

/n„n ;3 2a \ / Q 2 -L 2 ' 
-So/, - 2 - — c lW n a) n->) I M — 



r J r 2 \ 2r 



! )\ z A 



Yr — 



! )\ z 2 I : 

1 ■ — 



- ■ / \ / ' J '' 

2a x / Q 2 -L 2 \ . a 2 
-2^3V>n a) rF) yM — 

( Q 2 -L 2 

+2 M - 

\ 2r 



s o ( n * n P 2a 
i? a/3 = y— — 



+ 



Q 2 -L 2 



J e^e^rfn 1 - 4 ( M — — J -n^S a 

40(Q 2 -L 2 )\ z 4 \l f/ 



! )V 2 \ 



2 



r s : (- 2 + ^) 
'l4M^ 







; 2 / 40 (Q2 _ L : 
-+ I63M ^_ 



a 4 



/ 5(Q 2 -L 2 )\^ 2 / 4(Q 2 -^ 

8M - — + 7M - 

y r y y r 



12 (m-«1^!) 

[ \ 2r ) 

1) - (tm- 4 w 2 ; l2 >) 



7^ IV Y* } \ T J f 



—± \ 8 p<* 



(31) 
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and the asymptotic form of the electromagnetic potential has the form 

A °-^°{^)' ^=-£S +0 -(p)- (32) 

The determinant of the metric tensor is 



9 = ~l- 



(33) 



Now we evaluate all the required components of the contorsion and the basic vector using 
((2) and (4)) neglecting terms beyond the fourth powers of a. 

Substituting Eq. (31) and all components of the contorsion and basic vector in (24) we 

get 



2n c 



(M 



Q 2 -L 
2r 



2 a 2 



-A M ( l --^)-^——{i-^r, 



8r 2 



' f, z 2 189z\ 

M (9-126- + —— ) 



Q 2 - L 2 \ . 9Qz 2 120*V 



+ 



Az5° 



KT 



a 



(M — 



Q 2 -L 2 



21z 2 



2r 



+ in M(9 -— » 



Mnpe a/33 (a a 



(34) 



and 

/ 0/3 _ 
4a — 



2kt 2 




r 


2(g 2 - 






r 




35(Q 2 


-L 2 ) 



2r 



2 fi 5(Q 2 - L 2 ) 



zr 

(- - ^) «' + [(- - - h - 2J ^ 1 ) 3 

+ (lOM - V *) - (5M - *M) ^} + [( 2 M - 

V 8r J r 2 \ r ) r A \ 8r J r 6 

(1 15 ^!) W ) + + («- - *«£z£>) £! - ( 120 m 



e a3 /J + 2^3M 
5/ - (1 - 3^^)5 Q V) - (24M 
^ £a73 ^nV + 2 ( 8 M - 

' ! )\ : 2 ~ 



2 )\ z 



2r 

! )\ ( 2A(Q 2 -L 2 " 
' - ( 35M ^- ' 



2M - 

4r 



i a "-(i-i5^)« a vi + 

693(Q 2 -L 2 ) ^ - 4 " 



8r 



I — 

I r 4 



+ fl20M- 693 ' Q2 - 
\ 8r 



e^e^raV + 2 



/ 7*^ V 

W- 105 < g2 - L2 ))-f96M- 315 < g2 - ir 



J )\ ^ 2 



L 2 )\z 4 



^n% 3 + 



8r ; v 
12M _105(Q*-L*)V 
8r / r 2 



„, 693 Q 2 -L 2 \z 6 

120M ^ ' 

\ 8r 



(35) 
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The components we are interested in of the energy-momentum density are given by 



T 



o _ Q 2 - L 2 



KT 



a 2 f 6 . , ,.\ a 4 / 24, 9 9 , 
1 -^( 2 -^ 2 + ?/ 2 )J + 74( 3 -^ 2 + ^ 



= 2a(Q 2 - L 2 )e a/33 r/ / + 3a 2 _ 2 

Kr 5 ^ r 2 



30 



2 i 2\2 

4 + y ) 



(36) 



Further substituting (36) in (28) and then transforming it into spherical coordinate we 
obtain 



E = 



(Q 2 - L 2 ) 



KT 



l l l + 4 ° 2r2 ~~ 6fl2r2 C ° S2 9 + ^ ~ 36tt4 C ° S2 ^ + 30a4 C ° S4 B ) ^ Sil1 

(37) 



Performing the above integration we get the energy associated with the exterior general 
axisymmetric black hole given by Eq. (22) in the form 



E(R) 



total 



c.i terior {Q ^ ) / 1 , ^ \ 



R 



2 + 3i? 2 + 10i? 4 j + ° 



R 6 



(38) 



Eq. (38) shows that the total energy associated with the spacetime given by Eq. (22) is 
shared by its exterior as well as interior. 

Now we turn our attention to the angular momentum of the solution given by Eq. (22). 
From (36) we can get the components of the momentum density in the form 



KT 



(39) 



from which we can show that there is no momentum-density associated with Kerr black hole 
i.e,. when Q — L — 0. The components of the angular momentum of a general-relativistic 
system is given by [40] 

J n = f f l{x p P^ - x y Pfs)d 3 x , (40) 



where a, j3, 7 take cyclic values 1,2,3. Using (39) in (40) and transforming the expressions 
into spherical coordinates 



2a(Q 2 - L 2 



KT 



(r 2 + 3a 2 (1 - 2 cos 2 6) sin 3 OdrdOdcf) , (41) 



performing the above integration between two spheres of radii R\ and R2 we get 

R2 



J a = -2a(Q 2 - L 2 )e 12a 



3r 5r 3 



(42) 



Ri 



Eq. (42) shows that the rotation of the charged object is responsible for the angular mo- 
mentum distribution due to the electromagnetic field and the NUT parameter. 
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5. Main results and discussion 



In this paper we have studied the coupled equations of the gravitational and electromag- 
netic fields in the tetrad theory of gravitation. Applying the most general tetrad (15) to the 
field equations (7)~(9) we have obtained an exact solution given by Eq. (22). Eq. (22) is 
a general axi-symmetric solution from which we can generate the other known solutions like 
Schwarzschild, Reissner-Nordstrom, Kerr and Kerr-Newman spacetimes. The metric tensor 
associated with Eq. (22) is the Kerr-NUT spacetime. 

It was shown by M0ller [41] that the tetrad description of the gravitational field allows a 
more satisfactory treatment of the energy-momentum complex than does GR. Therefore, we 
have used the superpotential (24), formulated within the framework of teleparallel spacetime 
to calculate the energy- momentum density given by Eq. (27). We have used this energy- 
momentum density to evaluate the total exterior energy of the gravitating charged rotating 
body given by Eq. (22). Because the definition of energy of Eq. (28) requires its evaluation in 
Cartesian coordinate, the calculations without any approximation is obviously very tedious. 
Moreover, the intrinsic rotation parameter a is quantitatively very small for most physical 
situations. Therefore, for our convenience we keep terms containing powers of a up to the 
fourth order. With this approximation we have calculated the components of a covariant 
and contravariant tetrad fields (29) and (30), a covariant components of metric tensor (31). 
Calculating all the necessary components of the contorsion and basic vector ((2) and (4)) 
and using (31) in (28) we have obtained the expression of the exterior energy of the general 
black hole given by Eq. (22) till the fourth order. When the rotation, the charge as well as 
the NUT parameters are considered we get an additional terms 



which is the energy of the exterior magnetic field due to the rotation of the charged object. 
The asymptotic value of the total gravitational mass of a Kerr-NUT spacetime is the ADM 
[42] therefore, the energy associated with a Kerr-NUT contained in a sphere of radius R is 



Switching off the rotation and NUT parameters we found that the energy given by Eq. (44) 
will be the same as of Reissner-Nordstrom metric [43, 44]. The energy of Eq. (44) is confined 
to its interior only when we set the charge and the NUT parameters to be zero. This result is 
quit in conformity with those of Virbhadra [45, 46] and Cooperstock et al. [47] and Ahmed 



Using the value of energy-momentum density given by Eq. (36) we have calculated the 
momentum density. Also we have calculated the angular momentum distribution due to 
the electromagnetic field present in the Kerr-NUT field using Eq. (39) in Eq. (40). It is 
clear from Eq. (42) that the momentum density depends mainly on the rotation of the 
charged object and has only one component. As is clear from Eq. (42) that the angular 




(43) 




(44) 



[48]. 
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momentum depends on the even powers of the charge and NUT parameters and odd powers 
of the rotation parameter which means that the direction of the angular momentum vector 
depends on the direction of the rotation parameter. 
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